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Off-diagonal kinetic theory in ultrasoft momentum region at high temperature 
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We derive a linearized kinetic equation for fermionic excitations with an ultrasoft momentum, 
g 2 T, from the Kadanoff-Baym equation in a Yukawa model and quantum electrodynamics (QED) 
at extremely high T, where g is the coupling constant. We show that this equation is equivalent to 
the self-consistent equation in the resummed perturbation theory used in the analysis of the fermion 
spectrum with the ultrasoft momentum. Furthermore, we derive the equation that determines the 
n-point function with external lines for a pair of fermions and (n — 2) bosons with ultrasoft momenta 
in QED. 
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I. INTRODUCTION 

Rclativistic plasmas at extremely high temperature T 
such as electron and/or quark-gluon plasma are expected 
to be realized in the early universe and the experiments 
of heavy ion collisions at Rclativistic Heavy Ion Collider 
and Large Hadron Collider. These systems have multi- 
energy scale at weak coupling (j < 1) even at so high 
T that all mass scale in the Lagrangian is negligible; e.g, 
typical energy of particles is of order T, while the en- 
ergy of the bosonic and fermionic collective excitations 
called plasmon [1] and plasmino [2] is of order gT, and 
the fermion damping rate in gauge theories is of order 
g 2 T [3]. These energy scales, T, gT, and g 2 T, are called 
hard, soft, and ultrasoft scales, respectively. 

Unlike perturbation theory in the vacuum, the loop 
expansion at finite temperature does not correspond to 
the coupling nor h expansions, and the scale dependent 
expansion is necessary (Fig. 1). When the energy scale 
is gT, the one-loop approximation obtained by integrat- 
ing out hard internal momenta, called hard thermal loop 
(HTL) approximation [4], is reliable, and the approxi- 
mation establishes the HTL effective theory [5]. If the 
energy is of order or much less than g 2 T, some resum- 
mation 1 is necessary [6-12]. Here let us focus on the 
analysis for the fermion propagator with an ultrasoft mo- 
mentum. The rcsummation scheme used in the analysis 
of the fermion propagator consist of the following proce- 
dures [6, 8]: 

1. Resumming the thermal mass and the decay width 
of the hard particles. 
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1 This rcsummation [6-12] is different from the well-known HTL 
rcsummation, which is the resummation for the contribution 
from soft-internal momenta [5]. 



2. Summing all the ladder diagrams (in gauge theory). 

By analyzing the fermion propagator using the resummed 
perturbation theory, a novel fermionic excitation was 
found [6, 8]. 

Each perturbation scheme with soft and ultrasoft mo- 
menta can be interpreted as different kinetic equations. 
When the energy is soft, the effect of collisions is negli- 
gible and hence the collisionlcss kinetic equation called 
Vlasov equation is valid [13]. This equation corresponds 
to the HTL approximation [4] introduced above. When 
the energy is of order or less than g 2 T, the effects of the 
interaction among the hard particles are no longer neg- 
ligible in general, so it is necessary to take into account 
their effects. In fact, the analysis on the gluon self-energy 
with an ultrasoft momentum should be performed with 
the resummed perturbation theory including the effect of 
collisions, and it was done by constructing the linearized 
Boltzmann equation [7]. The computation of transport 
coefficients, whose energy scale is much smaller than g 2 T, 
is also performed with the resummed perturbation the- 
ory [10], the linearized Boltzmann equation [11], and the 
two-particle irreducible formalism [12]. These methods 
produce the same result in the leading order of the cou- 
pling constant. Using the correspondence between the 
perturbation theory and the linearized kinetic theory, 
the self-consistent equation in the resummed perturba- 
tion theory is interpreted with the language of the kinetic 
theory. 

In this paper, we derive a generalized and linearized 
kinetic equation for fermionic excitations with an ultra- 
soft momentum, which wc will call "off-diagonal" kinetic 
equation later, in the Yukawa model and QED, while 
the Boltzmann equation discussed above is employed for 
bosonic excitations. Our equation is systematically de- 
rived from the Kadanoff-Baym equation, and is equiv- 
alent to the self-consistent equation in the resummed 
perturbation theory [6, 8] used in the analysis of the 
fermion propagator. The derivation helps us to establish 
the foundation of the resummed perturbation scheme. 
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The kinetic equation will also give us the kinetic inter- 
pretation of the resummation scheme. Furthermore, we 
also discuss the procedure of analyzing the higher point 
functions not only two-point function of the fermion in 
QED. 

This paper is organized as follows: Section II is devoted 
to the derivation of the generalized and linearized kinetic 
equation and the discussion on the kinetic interpretation 
of the self-consistent equation in the resummed pertur- 
bation theory in the Yukawa model, which is the simplest 
fermion-boson system. In Sec. Ill, a similar analyses in 
QED is done in the Coulomb gauge. We make conclud- 
ing remarks and summarize our results in Sec. IV. In 
Appendix A, we evaluate the structure of the induced 
source terms. We analyze the Kadanoff-Baym equation 
that would yield the linearized kinetic equation in the 
diagonal case in Appendix B. Appendix C is devoted to 
showing that the temporal component of K^{x,y) intro- 
duced in Sec. Ill can be neglected. We show that the 
result obtained in Sec. Ill is valid also in the temporal 
gauge, which is the gauge-fixing of the original paper [8] 
of the resummed perturbation theory, in Appendix D. 



II. KINETIC EQUATION IN YUKAWA MODEL 

In this section, we derive a novel linearized ki- 
netic equation from the Kadanoff-Baym equation in the 
Yukawa model, which is the simplest model for fermion- 
boson systems. We will find the vertex correction is neg- 
ligible, which makes the analysis simpler than that in 
gauge theories. Next, we show that the kinetic equa- 
tion is equivalent to the self-consistent equation in the 
resummed perturbation theory [6, 9], and discuss the in- 
terpretation of the self-consistent equation using the cor- 
respondence to the kinetic theory. 



A. Derivation of the kinetic equation 

Throughout this paper, we work in the closed-time- 
path formalism [14, 15]. We perform the derivation of 
the kinetic equation in a similar way used in [13, 14] by 
applying the gradient expansion to the Kadanoff-Baym 
equation [16] and taking into account the interaction ef- 
fect among the hard particles in the leading order. 

Let us consider the following situation to analyze the 
fcrmionic ultrasoft excitation: Before the initial time 
to, the system is at equilibrium with a temperature T. 
Then, a (anti-) fcrmionic external source rj(x) {rj(x)) and 
a scalar external source j{x) are switched on. As a result, 
the system becomes nonequilibrium. We will consider the 
case that j(x) and rj(x) vanish and r){x) is so weak that 
the system is very close to the equilibrium, i.e., the lin- 
ear response regime. Concretely, we will retain only the 
terms in the linear order of the fermionic average field W 
in the fermionic induced source, which will be introduced 
later. 



Let us consider the generating functional in the closed 
time formalism [14], 



(2.1) 



with 



S= / d 4 4£[0,^]-(j0 + V^+^)], (2.2) 
Jc 

where <p and tp are the scalar and the fermion fields. The 
space-time integral is defined as J c d A x = J c dx° J <i 3 x, 
where J c dx° is the complex-time integral along the con- 
tour C = C+ U C~ U C° in Fig. 2. We will take t -> -oo 
and tf — ^ oo to factorize out the contribution from the 
path C . The Lagrangian in the Yukawa model with the 
massless fermion and boson has the form, 



(2.3) 



where we did not take into account the self-interaction of 
</> such as A^> 4 for simplicity. By performing an infinites- 
imal variation with respect to <j) or ip in Eq. (2.1), we 
obtain the following equations of motion: 

ip x [$]¥(x) = v(x) + Tfind(x), (2.4) 
-d 2 $(x) - gW{x) = j(x)+j ilLd (x), (2.5) 

where # = ((/>) (<& = (ip) ) is the expectation value of the 
scalar (fermion) field, and -pz^] = <j) x + ig^{x). Here 
the expectation value for an operator O is defined as 



(O) = ^ / V(t>V%pV^e tS 0. 



(2.6) 



mnd(x) = g(4>(x)ip(x)) c (Jiad(x) = g{ip{x)ip{x)) c ) is the 
fermionic (scalar) induced source, and the subscript c 
denotes "connected," i.e., (cf>(x)ip(x)) c = (cf>(x)ip(x)) ~ 
$(x)\P(x). 

By differentiating Eq. (2.4) with respect to j(y) and 
Eq. (2.5) with respect to rj(y), we obtain 

ip x [$]K(x, y) - gD(x, y)*{x) = (2.7) 



d 2 x K(y, x) - g(W(x)(^(yWx)) c + S(y, x)V{x)) 

_ . 5jind{x) 

$rj(y) 



(2.8) 



Here we have introduced the following propagators: 

(2.9) 
(2.10) 



D(x,y) = (T C (/>(x)</)(y)) c = *f^y> 
S(x,y)^(T c ^(xmy)) c = t S S y) 



Srj(x) ' 

K(x, y) = (Tc^(x)^(y)) c = = t -*g!() (2 .n) 

5j{y) 6r]{x) 
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FIG. 1. The correspondence between the resumed perturbation theory and the kinetic theory, and the fermionic modes obtained 
from these theories. The vertical axis denote the energy scale. 
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FIG. 2. The contour path C in the complex x° plane. 



where Tc means the path ordering on the complex-time 
path C; explicitly, 

D(x,y) = ec(x ,y a )D > (x,y)+6 c (y ,x )D < (x,y), 

S(x, y) = 6 c (x°, y Q )S>(x, y) - 6 c (y°, x°)S<(x, y), 

K(x, y) = 6 c (x ,y )K>(x, y) + 9 c (y°, x°)K<(x, y), 

(2.12) 

with 



D>(x,y) 


= {<j>{x)(j>{y)) c , 


(2.13) 


D<(x,y) 


= (cf)(y)(j){x)) c , 


(2.14) 


S>(x,y) 


= (ip(x)ip(y)) c , 


(2.15) 


S<(x,y) 


= (i/)(y)ip(x)) c , 


(2.16) 


K>(x,y) 


= (t/j(x)<p(y)) c , 


(2.17) 


K<(x,y) 


= (<Ky)ii>(x)) c , 


(2.18) 



and 6c(x,y) being the step-function along the path C. 
In the approximations introduced later, we can see that 



K > {x^y) and K < (x, y) coincide, which can be checked by 
K R (x,y) = i9(x°,y°)[K>(x,y)-K < (x,y)] ~ 0. For this 
reason, we simply write these two functions as K{x, y) 
from now on. We call K(x, y) "off- diagonal propagator" 
which mixes the fermion and boson, while we call -D(x, y) 
and S(x, y) "diagonal propagators." As will be seen in 
Sec. IIC, in the calculation of the ultrasoft fermion self- 
energy, the off-diagonal propagator is more relevant than 
other diagonal ones. 

By setting x° G C+ and y° e C~ in Eqs. (2.7) and 
(2.8), we obtain 



ip x [$]K(x,y) - gD<(x,y)V(a 



.^indO) 



dpv(x, y )+g(W(y)(^(y)iP(x)) c + S < (x,yp(y)) 

■Sjind(y) 



(2.19) 



5r](x) 



(2.20) 



Here we have interchanged x and y in the second equa- 
tion. Let us evaluate the right-hand side of Eqs. (2.19) 
and (2.20) using the chain rule: 



5r] ind (x) 5&(z) Sr] ind (x) 5$(z)\ 



c d Z { 8*{z) Sj(y) + 5<P{z) 6j(y)) 
d A z{S{ Xl z)K{z, y) + S(x, z)D(z, y)), (2.21) 



<c 
drj(x) 



d 4 z l 6$(z)5j ind (y) t 8&(z) 5j ind (y) 



K 5rj(x) 5$(z) Srj(x) 5<P(z) 
d 4 z(n(y, z)K(x, z) + S(x, z)3{z, y)). (2.22) 

c 

Here we have dropped (Srji nd /5^)(5^/Sj) and 
(<5 j i nd (S'F /Sff) since they contain more 
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than one fermionic average field. We have 

also used E(x,y) = —iSr]i nd (x)/S'I r (y) and 
II(x,y) = -i8j ind (x)/5<P{y) = 77(y,x), where 
£ (77) is the fcrmion (scalar) self-energy [7, 14]. 
We also introduced the off-diagonal self-energy, 
E(x,y) = -i8i lind {x)/8<P{y) = -i5j ind {y)/5W{x). 

The self-energies are decomposed for arbitrary a; and 
y° on the time path C: 

n(x, y) = 9 c (x ,y°)n>(x, y) + 9 c (y°, x°)II<(x, y), 

(2.23) 

S(x, y) = 6 c (x ,y°)S>(x, y) - 9 c (y°, x°)S<(x, y), 



(2.24) 

E(x,y) = 6c(x O ,y O )S>(x,y) + c (y o ,x o )S<(x,y). 

(2.25) 

We have not taken into account contact terms, which is 
negligible in the leading order as will be seen later. From 
Eqs. (Al), (A2), (A7), and (A8) in the Appendix, we 
can rewrite Eqs. (2.19) and (2.20) in terms of real time 
integral instead of that on the complex-time-path: 

ip x [$]K(x,y) - gD < (x,y)f(x) 

/OO 
d A z(Z R (x, z)K(z, y) + S R (x, z)D<(z, y)), 
-OO 

(2.26) 

-d 2 K(x,y)+gS < (x,y)V(y) 

d 4 z(n A (z, y)K(x, z) - S<(x, z)H R (z, y)), 

(2.27) 

where the retarded fermion (advanced scalar) self-energy 
E R (x,y) = i6{x Q ,y°)[Z>{x,y) + E<(x,y)} (II A (x,y) = 
~i6(y a ,x°)[II > (x,y) — II < (x,y)]) and the retarded off- 
diagonal self-energy E R (x,z) = i8(x° , z°)[E > (x, z) — 
E < {x 1 z)] have been introduced, and the limits to 
— oo and tf — > oo have been taken. Here we dropped 
'I'(y)('ip(y)'ip(x)} c because {ip{y)ip{x)) c contains more 
than one 9. Equations (2.26) and (2.27) are the 
Kadanoff-Baym equations from which the kinetic equa- 
tion is derived. 

Let us introduce the Wigner transformation to derive 
the kinetic equation, which is defined as 



f(k,X) 



d*se iks f (X + -,X 



(2.28) 



where s = x—y, X = (x+y)/2, and f(x, y) is an arbitrary 
function. After performing the Wigner transformation, 
Eqs. (2.26) and (2.27) become 



'>x 



ig<P{X) K(k, X) + igD < {k, X)V(X) 



= i(-S R (k 1 X)K(k,X) - S R \k,X)D < {k,X)), 

(2.29) 

(k 2 - ik ■ d x )K(k, X) + gS< (fc, X)&(X) 

= n A (k,X)K(k,X) - S<(k,X)S R (k,X). (2.30) 



Here we have used the following transformation law un- 
der the Wigner transformation, 



f(x)g(x,y)^ f(X)g(k,X) 

-\{dxf)-{d k g) + ..., (2.31) 

f(y)g(x,y)^ f(X)g(k,X) 

+ \{dxf)-{d k g) + ..., (2.32) 
d zg(x, z)h(z, y) — ¥ g(k, X)h(k, X) 



+ ^{g,h}p.B. + ..., (2.33) 

where {g, /i}p. b. = d k g ■ dxh — dxg • d k h is the Poisson 
bracket, and neglected higher-order terms that contain 
dx since we focus on the case that the inhomogeneity of 
the average field is dx ~ g 2 T, while a typical magnitude 
of k is of order T. This expansion is called gradient ex- 
pansion [7, 13, 14]. We retained the second terms in the 
left-hand sides of Eqs. (2.29) and (2.30) because the first 
terms, which seem to be the leading terms in the gradient 
expansion, will cancel out in the next manipulation. 

By multiplying Eq. (2.29) by (-ift + $ x /2 + ig@(X) + 
iE R (k, X)), adding Eq. (2.30), and setting <7> = 0, we get 



(2ik ■ d x - S R {k, X)} + n A (k, X))K(k, X) 
= g(fiD< (k, X) + S< (k, X))P(k, X). 



(2.34) 



Here we have introduced gt(k,X) = g^(X)+~ R {k,X). 
We have neglected higher order terms of order g 4 T 2 K 
and g 3 T~ 1 r. In the leading order, the coupling depen- 
dence in D < (k) and S < (k) is negligible, so that D < (k) 
and S < (k) are replaced by the propagators at equilibrium 
and free limit (g = 0): 

D"<(k) = p°(fc)n s (fc°), (2.35) 
S 0< (k) = ^p°(k)n F (k°), (2.36) 

where p°(k) is the free spectral function given by 



p°(k) = 27rsgn(fc°)<5(fc 2 ) 

= |g(^ -|k|)-^ 



|k|)), 



(2.37) 



and n B (k°) = (e k ° ' T - l)" 1 (n F (k°) = {e k ° ' T + l)" 1 ) is 
the boson (fermion) distribution function at equilibrium. 
We note that though the massless condition k 2 = ap- 
pears in Eq. (2.37) in the present approximation, k 2 is 
expected to be of order g 2 T 2 if one takes into account the 
interaction at equilibrium. For this reason, we will use 
the order estimate k 2 ~ g 2 T 2 . We also note that K can 
not be replaced by that at equilibrium since K vanishes 
at equilibrium. 

We see that k 2 terms in the left-hand side of Eq. (2.34) 
were canceled out and k-dx ~ g 2 T 2 term remains. Thus, 
we can neglect the terms much smaller than g 2 T 2 K in the 
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calculation of the leading order. Following this line, the 
diagonal self-energies are replaced by those at equilibrium 
in the leading order, whose diagrams are shown in Figs. 3 
and 4: 

{ft ) E R ^(k)}=m 2 f , (2.38) 
n A ^(k) =ml, (2.39) 

where m/ = gT/(2y/2) and nib = gT/^/6 are asymptotic 
thermal masses for the fermion and the scalar boson in 
the leading order [17]. Note that the imaginary parts of 
the self-energies ~ g 4 T 2 / g) and momentum depen- 
dence are negligible due to higher order of the coupling 
constant. We have used the on-shell condition, k 2 ~ 0, 
which will be verified later. 

The same line as the diagonal self-energies case justifies 
substituting the off-diagonal self-energy in the leading 
order, shown in Fig. 5 diagrammatically. The expression 
is given by 

~ R (k,X)=-g 2 J ^L S 0R (k + k')K(k',X). (2.40) 

Here S 0R (k) = — ft/((ko + ie) 2 — k 2 ) is the free fermion re- 
tarded propagator at equilibrium. We note that the self- 
energies can not be neglected unlike dx ~ gT case 2 [13], 
because {ft,S}K, IIK, T S is the same order as 
g 2 T 2 K. 

Using these expressions, Eq. (2.34) becomes 
(2ik-d x + Sm 2 )K(k,X) 

(2 41) 

= gftp°(k)(n B (k ) + n F {k°))f(k,X), 

where Sm 2 = m 2 — m 2 . We note that K(k,X) becomes 
finite only when k 2 = because of 8(k 2 ) in the right- 
hand side. We also note that fiK(k,X) ~ g 2 TK(k,X), 
which is confirmed by multiplying Eq. (2.41) by ft from 
the left. This property makes the vertex correction term, 
ftS R (k, X), negligible, which corresponds to the fact that 
there is no vertex correction in the analysis using the 
resummcd perturbation theory [6, 9]. Thus we get 

(2ik-d x +6m 2 )K(k,X) 

= gftp°(k)(n B (k ) + n F (k a ))V(X). 

The schematic figure of K(k,X) is depicted in Fig. 6. 
The solid (dashed) line with the blob stands for the re- 
summcd fermion (boson) propagator. 

B. Kinetic interpretation 

By introducing the 11 off- diagonal density matrix" 
A±{k,X) defined &sK(k,X) = 2Tr5{k 2 )(0(k o )A + (k,X) + 



2 This is because k ■ d x K ~ gT 2 K > {ft, E}K, IIK, T^S ~ 
g 2 T 2 K. 




FIG. 3. The fermion retarded self-energy E R ^(k) in the 
leading order. The solid line is the fermion propagator and 
the dashed line is the scalar (the case of the Yukawa model) 
or photon propagator (the case of QED). We note that the 
photon propagator in the case of QED should be replaced 
with the HTL-resummed one [5] in this figure. 




FIG. 4. The boson advanced self-energy 77 A(cq) (fc) in the 
leading order. The notations are the same as Fig. 3. 



6(— fc°)T_ (— k, X)), we arrive at the following general- 
ized and linearized kinetic equation from Eq. (2.42): 

/ Sm 2 \ 

{^■^±^)a ± (Kx) (243) 

= gt(n B (\k\) + n F (\kW(X), 

where = (l,k) and k = k/|k|. We have to note 
that this equation is not a usual kinetic equation because 
A(k,X) can not be interpreted as a distribution func- 
tion since it is the propagator between different particles 
in the fermionic background. Nevertheless, we call this 
equation "generalized kinetic equation" because we can 
obtain the Boltzmann equation if we analyze the time- 
evolution of the diagonal propagator instead of the off- 
diagonal one [7]. In fact, Eq. (2.43) has the following 
points similar to the Boltzmann equation: 

• The particle is on-shell (k 2 ~ 0). 

• The equation has the common structure as 
the Boltzmann equation: both have the non- 
interacting part, the interaction part between the 
hard particle and the average ultrasoft field, and 
the interaction part among the hard particles, 
which correspond to the drift term, the force term, 
and the collision term in the Boltzmann equation, 
respectively. 

When dx ~ gT, Sm 2 is negligible, and Eq. (2.43) be- 
comes the counterpart of the Vlasov equation [13]. Let 
us recapitulate the interpretations of each term in Eq. 
(2.43) except for the dm 2 term. The first term in the 
left-hand side describes the time-evolution of A±(k 7 X) 
in the free limit {g = 0), so this term corresponds to 
the drift term in the Boltzmann equation. On the other 
hand, the term in the right-hand side expresses the effect 
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FIG. 5. The off-diagonal self-energy S(k, X) in the leading 
order. The propagator that is composed of the solid line and 
the dashed line with the black blob is the off-diagonal propa- 
gator. The other notations are the same as Fig. 3. 



O 



y 



FIG. 6. The schematic figure of the off-diagonal propagator 
K(x, y) in the Yukawa model in the leading order in the linear 
response regime. The solid (dashed) line with black blob is 
the resummed fermion (boson) propagator that contains the 
information on the fermion (boson) self-energy £ (II). The 
gray blob represents the fermionic average field Of. 



from the average fermionic field. Hence this term corre- 
sponds to the force term in the Boltzmann equation. 

Sm 2 term has no counterpart in the usual Boltzmann 
equation, which describes the time-evolution of the di- 
agonal propagators, S(k,X) and D(k,X). Therefore we 
cannot obtain kinetic interpretation of that term in the 
usual sense. 

Now let us discuss the origin of the Sm 2 term. The 
origin of the term is -{ft,S R (k,X)} + II A (k,X) in 
Eq. (2.34). The real parts of those terms are of order 
g 2 T 2 while the imaginary parts are g^T In <? _1 , so the 
contribution in the leading order comes from the real 
parts. The difference of the real parts of the diagonal 
self-energies expresses the difference of the dispersion re- 
lations of the scalar and the fermion, so we call the dm 2 
term "mass difference term." We note that the "mass" 
here is not the bare one but dynamically generated one 
thorough the interaction among the hard particles. 

Here let us see the reason why the mass difference 
term docs not have its counterpart in the diagonal 
case. In the case of the kinetic equation on S(k,X), 
the corresponding term to the mass difference becomes 
-{#, S R (k, X)} + {ft, S A (k, X)} = -2ilm{#, S R {k, X)} 
instead of -{jk, E R (k,X)} + n A (k,X), so that the real 
part is canceled out (See Appendix B for the detailed 
derivation) . 

We note that this term, which is pure imaginary as 
a result of the cancellation of the real part, is of order 
<? 4 Tln(l/<7) and thus negligible since dx ~ g 2 T. We 
also see that the terms coming from the self-energy have 
complicated form in the diagonal case, while they are re- 
duced to the simple form, — {ft, E R (k)} + II A (k) term 
in Eq. (2.34), in the off-diagonal case in the linear re- 
sponse region. This difference comes from the following 




FIG. 7. The fermionic induced source ??i n d(^) in the leading 
order in the Yukawa model. This diagram is obtained by 
connecting two ends in Fig. 6. By truncating ty, we obtain the 
diagram expressing £ in the resummed one- loop analysis [6]. 



two facts: 

1. Since we linearize the equation in terms of the devi- 
ation from the equilibrium state, the terms contain- 
ing 5E and S° in Eq. (Bll) do not have its counter 
parts in the off-diagonal case [K(k,X) vanishes at 
equilibrium.] 

2. The structure of the right-hand sides of Eqs. (2.29) 
and (2.30) after neglecting the vertex correction 
terms are simpler than those of Eqs. (B4) and (B5). 
It is because K < (x 1 y) = K > (x,y). 



C. Correspondence between kinetic theory and 
resummed perturbation theory 

Here, let us show the relation between our kinetic equa- 
tion, Eq. (2.42), and the self-consistent equation in the 
resumed perturbation theory in Ref. [6]. For this pur- 
pose, we evaluate the retarded self-energy of the fermion, 
S R (p), which can be obtained from the linear response 
theory [13, 14] in momentum space as 



r, ind (p) = E r (pMp)- 



(2.44) 



On the other hand, ?7i n d can be written as ?7ind(^0 = 
gK(x,x) = gf d 4 k/(2n) 4 K(k 7 X) in our theory. Thus, 
the induced source is obtained from Eq. (2.42) as 



Vind(X) = g 2 



d A k ftp°(k)(n B (k°) + n F (k )) 



(2tt)' 



(2ik ■ d x + Sm 2 ) 



By performing the Fourier transformation, 



f(k,p)= I d*Xe^ x f{k,X) 



we obtain 



Vindip) = 9 1 



d 4 k ftp°(k)(n B (k°) + n F (k )) 



(2k-p + 6m 2 ) 



(2.45) 



(2.46) 



Hp)- 

(2.47) 
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Comparing Eq. (2.44) with Eq. (2.47), we obtain the self- 
energy, 



with 



d 4 k ltp (k)(n B (k°) + n F (k )) 



(2tt) 4 {2k-p + Sm 2 ) 



(2.48) 



This expression coincides with that of the retarded 
fermion self-energy in Ref. [6] except for the absence of 
the decay widths of the hard particles in the denomina- 
tor. As mentioned in the previous subsection, the decay 
widths of order g 4 T\n(l/g) is neglected when the exter- 
nal momentum is of order g 2 T; one can include them by 
taken into account the imaginary part of Eq. (2.42) if one 
is interested in the decay width. 

The diagrammatic representation of the fermion re- 
tarded self-energy in our formalism is the same as that in 
the resummcd perturbation theory [6] , which is explained 
as follows: The off-diagonal density matrix yl±(k, AT), 
which follows the generalized kinetic equation, is repre- 
sented by Fig. 6. rji n d(X), shown in Fig. 7, is diagram- 
matically obtained by connecting the ends of fermion and 
boson propagators in the right-hand side of Fig. 6. This 
diagram is the resummed one-loop diagram appeared 
in [6] itself except for the fermion average field . 



III. KINETIC EQUATION IN QED 

First we introduce the background field method, which 
is useful to construct the equations for the average fields 
and the Kadanoff-Baym equation in a gauge-covariant 
form. Next, we derive the generalized kinetic equation in 
the linear response regime adopting the Coulomb gauge 
fixing, in which the transversality of the free photon 
propagator simplifies the analysis. After the derivation, 
we show the equivalence between the generalized kinetic 
equation and the resummed perturbation theory [6, 8], 
and discuss the interpretation of the terms in the kinetic 
equation. We also check that the Ward-Takahashi iden- 
tity for the self-energy, consequence of the U(l) gauge 
symmetry, is satisfied in our formalism. Finally we dis- 
cuss how to compute the higher-point-vertex function 
whose external momenta are all ultrasoft, and make an 
order estimate of it at the weak coupling regime. 



A. Background field gauge method 

In the derivation of the average field equation and the 
Kadanoff-baym equation in QED, it is convenient to for- 
mulate them in a covariant form under gauge transfor- 
mations. For this purpose, we employ the background 
field gauge method [14, 18]. In this method, the follow- 
ing generating functional is employed: 



Z\j,r),ri;A,*,V]= [Va\[V^][V^]e 



js 



(3.1) 



S = I d 4 x 



C[A^ + a M , & + V, & + 4'} + C 



GF 



(3.2) 



where we dropped the ghost term, which is not coupled 
with the other fields. and are vector fields, \P 
and ip (ijj) are (anti-) spinor fields, and j M is the external 
current, respectively. The Lagrangian of QED has the 
form, 



f'H^H +iipp[a}ip, (3.3) 



d u a^ is the field strength and 



C[a,i(j,ip] = -i 
where F^[a] = d, 

D^la] = + iga^ is the covariant derivative. We have 
used g as a coupling constant of QED instead of the stan- 
dard notation e to make it clear that the same order 
counting appears as that in the Yukawa model. In the 
background field method, the fields in the Lagrangian are 
decomposed to the classical field, identified as the average 
fields later, and fluctuations in Eq. (3.1). The external 
sources are chosen to be coupled to a M , ip, and ip, but not 
to A^, and <P. We impose the following conditions: 



(o M ) = W = (V) = 0, 



(3.4) 



which implies that A^ and & (i') can be interpreted as 
the average parts of the photon and (anti) electron field, 
respectively, and — In Z coincides with the effective ac- 
tion [14, 18]. 

In the background gauge field method, the gauge-fixing 



term is chosen to be a functional of a p such as 



£gf — —A 



(G[a]f 



(3.5) 



where G[a] is the gauge-fixing function and A is the 
gauge-fixing parameter. Although a M is fixed by the 
gauge fixing term, the generating functional, Eq. (3.1), is 
invariant under the background field gauge transforma- 
tions defined by 

W(x) -> h(x)V(x), A^x) -> A^x) - ^h(x)d fi h\x), 

ip(x) ->• h(x)ij)(x), ip(x) -> ip(x)h'(x), a^{x) ->• a^x), 

T)(x) -> h(x)j](x), rj{x) -> fj{x)h\x), f{x) -> f(x), 

(3.6) 

where h(x) = cxp[i9(x)]. 

Since the fluctuations covariantly transform under 
Eq. (3.6), the propagators also covariantly transform as 

D^(x, y) = (T c a M (a;)a I/ (y)) c -> D fll/ (x, y), (3.7) 
S(x,y)^h(x)S(x,y)rf(y), (3.8) 
K^x, y) = (T c ^(x)a^y)) c -> h(x)K^x, y). (3.9) 
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Also, the Wigner transformed off-diagonal propagator 
in the leading order of g is covariant, which can be con- 
firmed by performing the gradient expansion [14]: 

K ll (k,X)^h(X)K li (k,X), (3.10) 

which implies that the Kadanoff-Baym equation covari- 
antly transforms with respect to the background gauge 
transformations as will be seen in the next subsection. 

We note that, apart from the covariance with respect 
to the background field gauge transformation, the gauge- 
fixing dependence, which will be confirmed only in the 
Coulomb gauge and temporal gauge in this paper, should 
be analyzed. 

B. Derivation of the kinetic equation 

We work in the Coulomb gauge-fixing condition be- 
cause this gauge fixing makes the analysis simple owing 
to the transversality of the free photon propagator. The 
gauge-fixing condition is G[a) = did 1 and A — > oo, which 
constrain the off-diagonal propagator as 

dj l K i (x,v)=0. (3.11) 

The equations of motion for the average fields are given 

by 

ip x [A)V{x) = n(x) + Vind{x), (3.12) 
8 2 A f "(x) - d»d v A u {x) - gV(x)'fV(x) 

= j"0«0 +J&d(*)- ( 3 - 13 ) 

Here we have used Eq. (3.4), and the induced fermionic 
source and the induced current are defined as 

Vind{x) = g(4(x)ip(x)) c = gj^K^x, x), (3.14) 
j£Jx) = g@(x)^(x)) c = gTr(^S<(x,x)), (3.15) 

which transform 77i n d(x) — > h(x)r/ ln d{x) and j^{x) — > 
Jind( a; ) un der the background gauge transformations. 
Therefore, Eqs. (3.12) and (3.13) transform covariantly 
with respect to the background gauge transformation. 
The equations for the propagators are given by 

p x K»(x, y) + ig^D^ix, y)V(x) 

/OO 
d 4 z(i; R (x,z)Ki 1 (z,y) (3.16) 
-oo 

+ S R (x,z)D<^(z,y)), 
(d 2 g^ - d»d v ) y K v (x, y) + gS<(x, y)>f*{y) 

/OO 
d*z(n A ^(z, y)K„(x, z) - S<(x, z)~ R »(z, y)), 
-oo 

(3.17) 

where 77 ^(x, y) and S^x, y) are the photon and the off- 
diagonal self-energies, respectively. Here we set x° E C + 
and y° G C~ . 



The Wigner transformed equations read 

+ ig lv D<^(k,X)y{X) 

= ~i(Z R (k, X)K»(k, X) + 3*{k, X)D <v »(k, X)), 

(3.18) 

(-fc 2 + ik ■ d x )K^(k, X) + fc^ n 7\r (fc, X) 
+ gS<(k,X)r*(X) 

= n A ^{k,x)K,(k,x) - s < (k,x)s Rfl {k,x), 

(3.19) 

where Eq. (3.11), the estimation K° ~ g 2 K l , shown 
in Appendix C, and the gradient expansion have been 
used. Multiplying Eq. (3.18) by [-ift + $ x /2 + igfi(X) + 
i£ R (k,X)) and Eq. (3.19) by P^(k) defined below, and 
subtracting the latter from the former, we obtain 

{-2ik ■ d x + 2gk ■ A(X) + {ft, S R {k, X)})K l (k, X) 

+ pl(k)n A ^(k,x)K u (k,x) 

= -g(^D <m (k, X) + P Tm (k)S<(k, X))t(k, X), 

(3.20) 

where we have introduced the projection operator into 
the transverse component, P? v {k) = g^g^ji^ij — kikj) 
, and gf^(k,X) = gj^(X) + » R »(k,X), and used 
pT(k)Kf*(k,X) = -Ki(k,X) = K l {k,X). One can 
show that the background fields and the coupling depen- 
dences in the diagonal propagators are weak 3 , so that 
we can replace the electron and the photon propagator, 
which is given as follows, by that in the free limit at 
equilibrium: 

D°<{k) =p°(k)n B (k°)P^(k). (3.21) 

The diagonal self-energies at on-shell in the leading order 
are given by 

S R{cq \k)} = m\- 2i( e k°, (3.22) 
pT(k)II A W a »{k) = -m 2 P?»(k), (3.23) 

where m e = gT/2 and to 7 = gT/\/6 are the asymp- 
totic thermal masses of the electron and the photon in 
the leading order. £ e ~ g 2 T\n(l/g) is the damping rate 
of the hard electron [3], which is much larger than the 
case of the Yukawa model. Unlike in the case of Yukawa 
model, the damping rate of the electron £ e cannot be ne- 
glected because ( e k° ~ inf.. while- the photon damping 



3 It can be shown by evaluating SS and 8D^ V in the diagonal 
kinetic equation: 8S ~ gA^ from Eq. (B12). Similarly, SD^ U ~ 
5S. These order estimates imply that A M contribution from S 
and in the right-hand side of Eq. (3.20) is much smaller 

than that from the left-hand side of that equation, by using the 
order estimate ~ g^ 1 ^. 
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rate of order g 4 T\n(l/g) can be neglected [19]. We note 
that the longitudinal part of the photon self-energy does 
not contribute because the projection operator P^ t (k) is 
multiplied. 

The off-diagonal self-energy in the leading order has 
the following expression, which has a similar form to that 
in the Yukawa model: 



S R »(k,X) 



d 4 k' 

W) 



: ys 0R (k + k')^K v (k',x). 

(3.24) 



The off-diagonal self-energy in r^(k,X) has to be re- 
tained in the case of QED unlike the Yukawa model since 
there is no special suppression mechanism. 

By substituting these expressions in Eq. (3.20), we ob- 
tain 



(-2ik ■ D x - 2iC e k° - £m 2 )iT(fc, X) 
= -g^P Tm (k)p°(k)(n B (k a ) + n F (k )) 
d A k' k a lu + 1 a k' v 



x n/„W(X)+g 



k ■ k' 



K a (k',X) 

(3.25) 



Here we have used P^ a {k)P^ a {k) = —P^ v (k) and intro- 
duced Sm 2 = m 2 — TOg. In this gauge- fixing condition, 
it is apparent that only the transverse component of the 
thermal photon contributes to if 4 (fc,Y) because of the 
projection operator P Tv ' l {k) appearing in the right-hand 
side of Eq. (3.25). We note that this equation transforms 
covariantly with respect to the background gauge trans- 
formation from Eq. (3.10). 

From Eq. (3.25), we write K % (k,X) in terms of the 
off-diagonal self-energy for later use: 

K\k, X) = g ^y^)(n^) + n,(*°) )fv 



2ik ■ D x + 2iC e k° + dm 2 



(3.26) 



The diagrammatic representation of this equation is 
shown in Fig. 8. 



C. Kinetic interpretation 

Next, we derive the linearized kinetic equation. Mul- 
tiplying Eq. (3.25) by ji from the left, we obtain 

{-2ik ■ D x - 2iQk° - 5m 2 )I/L'(k,X) 
= -g2^p (k)(n B (k )+n F (k°))<P(X) 
+ 7 < fPj < (ft)/j (fc)(n fl (ft ) + n F (k )) (3-27) 

(Pk' k a Y + k w ~i a 



x 9 



(2tt) 4 k ■ k' 



-K a (k',X). 



We decompose A±(k,X) into positive and negative en- 
ergy terms as K»{k,X) = 27r5(/c 2 )[0(fc o )yl+(k, X) + 



6(— k°)A^i(— k, X)], so that we arrive at the kinetic equa- 
tion from Eq. (3.27): 



2iv ■ dx ± 



Sm 2 



2iC e U±(k,X) 



= 2^[n B (|k|)+n F (|k|)]<F(X) 
- 5 2 7 ^K(|k|)+n B (|k|)]Pf (v) 



E 

s=± 



d 3 k' 1 s|kh Q 7„ ± Ik'lw 



(2tt) 3 2|k'| 



^ ±|K| ^ sa (k',Y), 



Ikllk'lu-t;' 



(3.28) 



where we have introduced z/ M = (l,k'). There are two 
terms that do not appear in the Yukawa model analyzed 
in Sec. II. One is the last term in the right-hand side. Be- 
cause the bare vertex term, which contains ^{X) and is 
in the right-hand side, is interpreted as the counterpart of 
the force term in the diagonal case [13], that vertex cor- 
rection term acts like 11 the correction to the force term," 
at least in the linear response regime. Note that this term 
mixes the positive and negative energy modes, unlike the 
case of the Yukawa model. 

The other is the third term in the left-hand side. This 
term has a similar form to the collision term in the re- 
laxation time approximation of the diagonal case, i.e., 
pure imaginary constant (2i£ e ) times A±(k,X). For this 
reason, we call this term "collision term." We note that 
this term is negligible in the case of the Yukawa model 
as shown in the previous section. 

In the diagonal case [7], the collision term contains 
momentum integral for the diagonal density matrix 4 . In 
contrast, in the off-diagonal case, such term does not sur- 
vive in the linearized equation. As a result, the collision 
term in the off-diagonal kinetic equation has a similar 
form to that in the relaxation time approximation. 

We emphasize that the off-diagonal self-energy is not 
negligible in the off-diagonal kinetic equation, while neg- 
ligible in the diagonal one because we neglect the higher 
order terms in W. This fact makes the correction to the 
force term, which is absent in the diagonal case, appears 
in Eq. (3.28). 

As we discussed in Sec. II B, both of the usual Boltz- 
mann equation and our generalized and linearized kinetic 
equation are composed of the non-interacting part, the 
interaction part between the hard particle and the av- 
erage ultrasoft field, and the interaction part among the 
hard particles. Which part is the counterpart of the mass 
difference term, the collision term, and the correction to 
the force term? Because the mass difference and the col- 
lision term come from the self-energies at equilibrium, 
they correspond to the interaction part among the hard 
particle. The correction to the force term is a part of 
the interaction part between the hard particle and the 
average ultrasoft field. 



The integral comes from the terms containing &E in Eq. (B12) 
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FIG. 8. The diagrammatic representation of self-consistent equation for K IJ '(x,y) in the leading order. For simplicity, is 
not drawn. 



TABLE I. The correspondence between the resummed per- 
turbation theory and the generalized and linearized kinetic 
equation. 



Diagrammatic method 


kinetic equation 


thermal mass difference 
decay width 
ladder diagrams 


mass difference term 
collision term 
correction to force term 



D. Correspondence between kinetic theory and 
resummed perturbation theory 



Here, let us show the equivalence between Eq. (3.25) 
and the self-consistent equation in the resummed pertur- 
bation theory [6, 8]. To this end, we rewrite Eq. (3.24) 
in terms of the off-diagonal self-energy using Eq. (3.26): 



r fl (k,x) = j^(x)-g 



d 4 k' 



(27T) 4 

lt'P^(k')p (k')[n B (k ,0 )+n F (k' )} 



2ik> ■ d x + 2< e fc'° + Sm 2 



x f a (k',X). 



(3.29) 



Here we set A M = 0. By performing the Fourier transfor- 
mation, Eq. (2.46), we get 



r»(k,p) = Y l -g 



d*k' 



(27T) 4 

^(^V(fc')K(fc' ) + n F (k'°)]r a (k', P ) 



2k' ■ p + 2i£ e k'° + Si 



(3.30) 



where r^{k,p) = r l± {k,p)W(p). We note that r il (k,p) is 
the vertex function introduced in [6, 8] whose momenta 
are hard and ultrasoft. Equation (3.30) is none other 
than the integral equation appearing in [6, 8]. 

The retarded fermion self-energy is also written in 
terms of the vertex function: from Eq. (3.26), we arrive 



at 

S R (p) 
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d A k 5${k,p) 



(2tt) 4 6&(p) 
d*k 7i ftPp(k)p (k)[n B (k ) + n F (k°)} 



(2tt)< 

xr»(k, P ). 



2k ■ p + 2i( e k° + Sm 2 



(3.31) 



This expression equals to that of the fermion retarded 
self-energy in the resummed perturbation theory [6, 
8]. Thus we sec that Eq. (3.25), derived in the 
non-equilibrium state in a linear response regime from 
the Kadanoff-Baym equation, is equivalent to the self- 
consistent equation in the resummed perturbation the- 
ory, which is constructed in the thermal equilibrium 
state. 

Here we discuss the correspondence between each 
scheme in the resummed perturbation theory [6, 8] and 
each term in the kinetic equation. As in the Yukawa 
model, the resummation of the thermal mass difference 
in the resummed perturbation theory corresponds to the 
mass difference term in the kinetic equation. The de- 
cay width corresponds to the collision term since both 
of them contains the damping rate of the hard electron, 
£ e • The ladder summation in the resummed perturbation 
theory [6, 8] is caused by the correction to the force term 
in the kinetic equation. Thus, the ladder summation cor- 
responds to the correction to the force term. These in- 
terpretations are summarized in the Table. I. 



E. Higher point functions 

The fcrmionic induced source r]i n a generates the higher 
point function, not only the fermion self-energy. In this 
subsection, we derive the self-consistent equation deter- 
mining the n-point function whose external lines consist 
of two fermions (iP) and (n — 2) bosons (A^) with ul- 
trasoft external momenta, and make an order estimate 
of the quantity. For example in the case of n = 3, we 
obtain the correction to the bare three-point function, 



Equation (3.30) is analytically solved for the energy region |po + 
*Ce| «9 2 T [6]. 
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91ii^{p -q-r), from r/ ind [13, 14]: 
5^(p-q~r)gSr^p,-q,-r) 



S^(q)SA^r) 
d 4 k' 
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A=0 

5 2 



(3.32) 



-6fi(k,p). 



(2tt) 4 8^{q)5A^{r) 

Here we have expanded K fl (k,X) around = 0: 

K"(k,X) = K> 1 (k,X) A=0 + 5K fl (k,X) + O{A 2 &), 

(3.33) 

where K fl (k, X)a=o contains one ^ while 8K^{k : X) con- 
tains one and one A^ . 

SK^ l (k, X) can be obtained by the following way. Col- 
lecting terms that contain one A^ in Eq. (3.25), we obtain 

(-2ik ■ d x - 2i( e k° - 5m 2 )6K l {k,X) 

+ 2gk-A(X)K t (k,X) A =o 
= ~g 2 $P TlJ \k)p\k)[n B {k a ) + n F (k )} 
d 4 k' k a lu + j a k' u 



(3.34) 



(2tt)< 



k ■ k' 



-SK a (k',X) 



Since K^i^k, X)a=o is determined by setting A^ = in 
Eq. (3.25), this equation is closed and 5K >1 (k,X) can be 
determined. 

Let us estimate the order of 6K^(k,X). From 
Eq. (3.25) and K% =0 ~ g^T^, we find 5K» ~ 
g- 1 T- 1 A^K% =0 ~ g- 2 T- 4 &A a . Therefore, the vertex 
correction is estimated as <7#-T M ~ which is much 

larger than the bare vertex, gj^ ~ g. Similar order es- 
timate for the n-point function with n > 3 can be done 
with the same procedure; result, we find that the 
order of the n-point- vertex function is g 2 ~ n . 



F. Ward-Takahashi identity 

Let us see that the off-diagonal self-energy given in 
Eq. (3.24) satisfies the Ward-Takahashi (WT) identity. 
From Eq. (3.24), we get 

k^(k,X) 



d 4 k' 



-I{(k',X)=g7 lind (X). 



(3.35) 



Here we have used $K v (k',X) = 0, which can be con- 
firmed by multiplying Eq. (3.25) by f from the left. This 
equation generates the WT identity derived in Ref. [6] by 
setting A 1 = 0. The WT identity implies that the vertex 
correction is not negligible because the identity relates 
the vertex correction to the fermion self-energy, which is 



much larger than the inverse of the fermion propagator 
with an ultrasoft momentum. In the Yukawa model, the 
WT identity associated with gauge symmetries is absent, 
so the smallness of the vertex correction is not in contra- 
diction with any identity. 

Equation (3.35) can be derived from the conservation 
law of the induced current, 



(3.36) 



By differentiating Eq. (3.36) with respect to \P(y), we 
obtain 



(3.37) 



Here S c (x~y) is the delta function defined along the con- 
tour C . By multiplying this equation by J d 4 s cxp(ik ■ s) 
and taking only the leading-order terms, we find 



ld {X) = -ik fi E R »(-k,X). (3.38) 



This is nothing but Eq. (3.35). The detail of the deriva- 
tion of Eq. (3.38) is shown in the Appendix A. 



IV. SUMMARY AND CONCLUDING 
REMARKS 



We derived the novel generalized and linearized ki- 
netic equation with an ultrasoft momentum from the 
Kadanoff-Baym equation in the Yukawa model and QED. 
Our kinetic equation is equivalent to the self-consistent 
equation in the resummed perturbation theory [6, 8] used 
in the analysis of the fermion propagator. This deriva- 
tion helps us to establish the foundation of the resummed 
perturbation theory. We gave the kinetic interpretation 
of the resummation scheme (Table. I). Furthermore, we 
made an order estimate of the higher-point function with 
ultrasoft external momenta, and obtained the equation 
determining the vertex correction in QED. 

In our analysis, the difference of the masses and the 
vertex correction, which reflects the fact that we analyze 
a process that changes the type of particle, plays an im- 
portant role. This fact suggests that these quantities may 
play an important role in the analysis of other ultrasoft 
quantity such as the flavor diffusion constant in QCD, 
where the masses of u, d and s are different. It is natural 
to ask whether such kinetic equation is derived in QCD 
because the investigation of the QCD at high tempera- 
ture is quite interesting since it is relevant to analysis of 
quark-gluon plasma [20], which is realized in heavy ion 
collision experiments. We hope to report the analysis in 
QCD elsewhere [21]. 
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Appendix A: ANALYTIC STRUCTURE OF 
INDUCED TERMS 



In this Appendix, we derive the right-hand sides of 
Eqs. (2.26), (2.27), and (3.38). First let us derive the 
right-hand sides of Eqs. (2.26) and (2.27): 



d zS(x, z)K(z, y) 



<fzE > (x,z)K{z,y)- / (fzS < {x,z)K{z,y) 



t°-ip 



d zS (x, z)K(z,y) 



d 4 z(S>(x,z) + S < (x,z))K(z,y) 



t°-%p 



d 4 zE < (x,z)K(z,y) 



i / d 4 zS H (x,z)K(z,y). 



(Al) 



In the last line we have taken t° —> — oo. We used the 
fact that the term integrated on C° becomes negligible 
in this limit [14]. In the same way, we get 



d i zll(y,z)K(x,z)~-i / d 4 zII A (z,y)K(x,z). 



(A2) 



Next, we evaluate the off-diagonal self-energy term. 



The vertex correction term of Eq. (2.21) becomes 

/ d i z~(x,z)D(z,y) 
Jc 



d 4 z~ > (x,z)D < (z,y) 
y° 

+ / d 4 z~<(x,z)D<(z,y) 

rt°-i/3 

d 4 z~ < {x,z)D > (z,y) 



d A z(E>(x,z)-E<(x,z))D<(z,y) 



(A3) 



+ / d i zS<(x,z)(D<(z,y)-D>(z,y)) 
Jt° 

'-i/3 

d i zS<(x,z)D>(z,y) 

X) 

d i z{E R (x,z)D < (z,y) 

+ E<(x,z)D A (z,y)), 

where the advanced boson propagator is D A (z 1 y) = 
—i9(y° — z°) [D > (z 7 y) — D< (z, y)} . Here we stop and dis- 
cuss the structure of the off-diagonal self-energy in the 
leading order in the (k,X) space. 

The off-diagonal self-energy in the leading order and 
linear response regime is given by 



£ (a;, y) = g 2 S°(x, y)K{y, x), 



(A4) 



where S°(x, y) is the free fermion propagator at equilib- 
rium. Thus, the components of S are given by 

S^(x,y)=±g 2 S°^(x,y)K(y,x). (A5) 

By performing the Wigner transformation, we get 

~^(k,X)=±g 2 J 0^S^(k + k')K(k',X) 7 (A6) 

with S a> (k) = fkp°(k)U - n F (k )). We see that since 
K(k',X) contains S(k' 2 ) and we focus on the on-shell 
case k 2 ~ 0, which are confirmed from the analysis 
in Sec. II A, (k + k') 2 ~ 2k ■ k' ^ 0. For this rea- 
son, S°^(k + k') ~ 0, which implies S^(k,X) ~ 0, 
so the only nonzero function of the off-diagonal self- 
energy appearing at Wigner-transformed Eq. (A3) is 
E R (k,X) ~ S A (k,X), Therefore, we drop the second 
term in Eq. (A3) because that term becomes negligible 
after the Wigner transformation, and hence the equation 
becomes 

/ d i zH(x,z)D(z,y)c- l ^ d 4 z~ R (x,z)D<(z,y). 

JC J-oo 

(A7) 
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In the same way, we get 



(TzS(x, z)~{z, y)~i d 4 zS < {x 1 z)E R (z, y). (A8) 



The calculation in QED can be performed in the same 
way. 

Finally we derive Eq. (3.38). By multiplying Eq. (3.37) 
by J d 4 scxp(ik ■ s), we get 

-gj d 4 S e lks ( ^^ { % (x) 



d 4 se ik - s d;S^y,x). 



Here we have set x°, y° G C + and neglected the sub- 
leading terms. The first term in the left-hand side has 
the same order of magnitude as the hard fermion self- 
energy S(k) times &{X), so that term is negligible. Thus 
the left-hand side becomes 77i n d(A). The right-hand side 
becomes 

/0 />oo 
d 4 se lk - s d s ^S > ^(y,x) + / d 4 se lk - s d s ^S < ^(y,x) 
-oo JO 

d 4 se ik - s df l E > ^(y,x) 



+ d i se ik - s d s JS<^y,x)-E>' i (y,x)) 



-ik^B >fl (-k,X) - k^S A >"(-k,X). 



(A10) 



We see that the first term in the last line is negligible be- 
cause of the on-shell condition. Thus we obtain Eq. (3.38) 
if we remember that E R (k,X) ~ S A (k,X). 



Appendix B: KINETIC EQUATION IN 
DIAGONAL CASE 



We derive the equation that corresponds to Eq. (2.42) 
in the diagonal case in this Appendix. The equation gov- 
erning the propagator of the fermion is as follows: 

PMS(x, y) + ig(K(y, x))^°W(x) 

Mnd(x) (Bl) 



= S c ^(x-y)+i- 



5r)(y) 



Since the second term in the left-hand side contains two 
'P, we neglect that term. Let us calculate the right-hand 
side. We set a; e C + and y° € C~ . Since the ver- 
tex correction term, which contains more than one i', is 



negligible, we obtain 
%nd(a;) 



5rj{y) 



d zS(x, z)S(z, y) 

d 4 zS > (x, z)S < (z, y) 
d 4 zS < (x, z)S < (z, y) 

.0 

d 4 zS < (x,z)S > (z,y) 



r" 

/ d 4 z(S>(x,z) + S<(x,z))S<(z,y) 



+ d 4 zE<(x,z)(S < (z,y) + S > (z,y)). 



t° 



(B2) 



By taking the limit t — > — oo, we get 
6rj(y) 

/OO 
d 4 z(S R {x, z)S < (z, y) + S < (x, z)S A {z, yj). 
-OO 

(B3) 

Here we have introduced the advanced fermion propaga- 
tor, S A (x, y) = -i6(y°, x°)(S>(x, y) + S<(x, y)). 
By performing the Wigner transformation, we get 



-it+^ + ig(*(X)-i?±.(d x ${X)) 



S < (k,X) 



i{E H (k, X)S (Jfe, X) + S< (k, X)S A (k, X)) 



(B4) 



The following equation is derived from the conjugate of 
Eq. (B4) by using j S>{x,y)j = S<(y,x): 



S<(k,X) 



if + ^-igU{X) + i^.(d x ^(X)) 



i(S R (k, X)S<(k, X) + S< (k, X)S A (k, X)). 



(B5) 



Here we have introduced the retarded fermion propaga- 
tor, S R , and the advanced fermion self-energy, S A , which 
are defined as follows: 

S R {x, y) = i6{x\ y°)(S>(x, y) + S<(x, y)), (B6) 
E A (x,y) = -ie(y°,x°)(E>(x,y) + S<{x,y)). (B7) 

By multiplying Eq. (B4) (Eq. (B5)) by -ifk + $ x /2 + 
ig {$(X) - id k ■ d x ${X)/2) - iZ R {k, X) (if + $ x /2 - 
ig{$(X) +id k ■ d x ${X)/2) - iS A (k,X)) from the left 
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(right), we get 

(-k 2 - ik ■ d x + gW<P(X) - i$d k ■ (dx$(X)) 
+ $(X)$ X ) + {$, E R (k, X)})S<(k, X) 
= -k , Z < (k,X)S A (k,X), (B8) 
S<(k, X)(-k 2 + ik-tx+ g(2ft$(X) + icF k ■ (d x $(X))ft 
-i%${X)) + {%£ A {k,X)}) 
= -S R (k,X)S < (k,X)ft. (B9) 

By subtracting Eq. (B8) from Eq. (B9), we get 

(2ik ■ d x - {ft, S R (k,X)}) S<(k,X) 
+ S<(k,X){ft,S A (k,X)} 
+ 2g$(X)[S<{k,X),fk} 
+ ig(d x $(X)){fk,d k „S<(k,X)} 

-i 9 0(X){^,d x s<(k,x)} 

= ft£ < (k,X)S A {k,X) - S R (k,X)S < (k,X)ft. 

(BIO) 

Here we linearize this equation. By introducing 
SS(k,X) = S(k,X) - S°(k) and 5£(k,X) = S(k,X) - 
E^ eq \k), we arrive at the following equation: 

hik ■ d x - {$, S (cq)R (k) - S^ A (k)}\ <5S*<(fc, X) 

- {^SS R (k,X)}S 0< (k) + S 0< (k){k',5£ A (k,X)} 
+ ig(d x <l>(X)){k',d ku S a< (k)} 

= k'5£ < (k,X)S 0A (k) - S 0R (k)SS < (k,X)^ 

+ ftE^ < (k)SS A (k,X) - 8S R (k,X)S ( - c ^ < (k)lt. 

(Bll) 

Similarly, we obtain the following equation in QED: 

(2ik ■ d x - {ft, S (cq)R {k) - E^ A (k)})5S < {k, X) 

- {ft,6E R (k,X)}S 0< (k) + S 0< {k){ft,SS A (k,X)} 
+ 2igk»d k -d x A li (X)S 0< (k) 
-^(X)[a^,S°<(k)] 

= k'5£ < (k,X)S 0A (k) - S 0R {k)SS < (k,X)f: 

+ ft£( cq '> < {k)5S A {k 7 X) - 6S R {k,X)U ( - c ^ < (k)l^. 

(B12) 



Appendix C: SMALLNESS OF K° 

Here we show that K° is negligible compared with the 
spatial components in the Coulomb gauge. For simplic- 
ity, we set = 0. We get the following equation by mul- 
tiplying Eq. (3.18) by -i$ + $ x /2 + iS R (k, X) from the 
left, subtracting Eq. (3.19) from the quantity obtained 



above, and setting fi = 0: 

(2ik ■ d x - {ft,Z R (k,X)} + (fc°) 2 )AVfc,X) 

-n Mv (k,x)K v {k,x) = g s < (k)r°(k,X) 

(fc°) 2 A' (fc, X) - n AOi (k, X)Ki(k, X) 

= gS<(k)f°(k,X). 

Since k° ~ T, we see that 

K° ~ g 2 K\ 

Appendix D: TEMPORAL GAUGE 



(CI) 



(C2) 



We show that the equation determining K M in the tem- 
poral gauge is the same as that in the Coulomb gauge, 
Eq. (3.25), in this Appendix because the resummed per- 
turbation theory was first proposed in the temporal 
gauge [8]. The gauge-fixing condition in the temporal 
gauge is G[a] = a and A — > oo. This condition is equiva- 
lent to the constraint a = 0. Because of this constraint, 
we have 

K (x,y)=D »(x,y)=0. (Dl) 
The equations governing K' are 

^+|+i 9 4(I)j K\k,X) 

+ ig^ j D <3i (k,X)&(X) (D2) 
= -i{E R (k,X)K i (k,X) + 5 R {k,X)D < i i {k,X)), 
(-k 2 + ik-d x )K l (k,X) 

*- l ^)\k-^yK j{k ,X) (D3) 
+ gS < (k,X) 1 l 1 , (X) 

= n M] (k, x)Kj (k, x) - s< (k, x)~ m (k, x). 

From these equations, we obtain 

(-2ik ■ d x + 2gk ■ A(X) + {ft, E R (k, X)})K\k, X) 

k^)\k-^) K A k,X) 

+ n AlJ (k 7 x)K J (k,x) 

= -{k'D <jl {k,X)+8 lJ S < {k,X))r 3 {k,X). 

(D4) 

Here let us evaluate k l Ki(k, X), which is the longitu- 
dinal component of K l (k,X). By multiplying Eq. (D4) 
by fci, we get 

|k| 2 (k- Kj(k,X) + k i n Aij (k,X)K j (k,X) 



-s < (k,x)k l r l (k,x). 



(D5) 
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Here we have neglected the terms that are of order 
g^^Kii^X). We see that WK^k^X) ~ g 2 TKi(k,X) 
and thus the longitudinal component of Ki(k,X), 
k z Ki(k,X), is negligible compared with the transverse 
component of Ki(k,X). We note that K° = 0, which 



is the result of the gauge-fixing condition, and k l Ki = 
are valid also in the Coulomb gauge in the leading order. 
Furthermore, also the free photon propagator at equilib- 
rium is the same as that in the Coulomb gauge. Thus, we 
can obtain Eq. (3.25) in the same way as in the Sec. IIIB. 
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